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Abstract 

Under hard-agar and nutrient-rich conditions, a cell of Bacillus subtilis grows as a 
single filament owing to the failure of cell separation after each growth and division 
cycle. The self-elongating filament of cells shows sequential folding processes, and 
multifold structures extend over an agar plate. We report that the growth process from 
the exponential phase to the stationary phase is well described by the time evolution of 
fractal dimensions of the filament configuration. We propose a method of characterizing 
filament configurations using a set of lengths of multifold parts of a filament. Systems 
of differential equations are introduced to describe the folding processes that create 
multifold structures in the early stage of the growth process. We show that the fitting 
of experimental data to the solutions of equations is excellent, and the parameters 
involved in our model systems are determined. 


1 Introduction 

Patterns observed in bacterial colonies growing on surfaces of semisolid agar plates realize 
a variety of fractal and self-affine structures studied in statistical mechanics and fractal 
physics [Tj. In a series of experimental studies [21 El H [5], it has been clarified that the 
morphology of growing bacterial colonies at the macroscopic scale does not depend on the 
biological details of individual organisms, but depends only on environmental conditions 
controlled by the agar concentration C a and the nutrient concentration C n . Cell motility is 
changed by varying C a and the growth rate is controlled by varying C n , and then different 
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patterns appear in different regions in a morphology diagram drawn on the C a -C n plane. On 
this plane, a diffusion-limited aggregation (DLA)-like pattern [BJ ;7j, an Eden-like pattern 
a concentric-ring pattern, a homogeneously spreading disk like pattern, and a dense 
branching morphology (DBM) pattern [T] have been recorded. Since some of these patterns 
are observed not only in biological systems, but also in chemical and physical systems, such as 
those in crystal growth, aggregation processes, and viscous fingering [10], the mechanism of 
pattern formation could be common in organisms and inorganic substances, and a theoretical 
study by mathematical modeling and analysis will be useful for understanding the underlying 
universal principles pn m m m, ng. 

In a recent paper [16], we have reported the physical aspects of the collective motion of 
bacterial cells observed in shallow circular pools prepared on the surface of an agar plate. 
The diameters of the pools are arranged so as not to be much larger than the length of the 
bacterial cells swimming in the pools. We used Bacillus ( B .) subtilis and found six different 
types of collective motion, including one-way and two-way rotational motion along the brim 
of a circular pool and collective oscillatory motion in the entire pool. Analyzing experimental 
observations in 117 circular pools, we found that these six types of collective motion can be 
classified using only two parameters: the reduced cell length A, which is defined as the ratio 
of the average cell length in a pool to the pool diameter, and the cell density p in the pool. 
We obtained a phase diagram for the collective motion drawn on the A -p plane and predicted 
that simple modeling with the two control parameters will be able to explain the variety of 
collective motion of bacterial cells. 

The above results imply that physical considerations are applicable and useful to explain 
the changes in the morphology of bacterial colonies at the macroscopic scale as well as 
the dynamical transitions of the collective motion of bacterial cells at the microscopic scale 
caused by environmental variations. On the basis of them, in this paper, we will report 
the experimental results and numerical analyses of growth processes starting from a single 
bacterial cell observed in formation of colony patterns. We tried to analyze the observed 
growth process by fractal analysis and by using systems of differential equations. 

Throughout the experiment reported in this paper, we used B. subtilis wild-type strain 
OG-01. Cells of this strain are rod-shaped (0.5-1.0 pin in diameter, 2-5 pm in length) with 
peritrichous flagella. They swim in a straight line in water by bundling and rotating the 
flagella. Under unfavorable environmental conditions, such as on a nutrient-poor medium 
or a dry agar plate, they become spores. When a small number of cells are inoculated on 
the surface of a semisolid agar plate, they go through a resting period of about 7 h before 
starting two-dimensional colony expansion. 

As briefly reported in previous papers , we observed string like objects in microscopy 
observations of an inoculation spot of a bacterial suspension in the later stage of the resting 
phase. In this study, we focus on the growth process that started from a single bacterial 
cell in the resting period under hard-agar and nutrient-rich conditions. In this case, cell 
multiplications are repeated with a constant cell cycle (doubling time), but daughter cells 
fail to separate after each cell cycle, although the cytoplasm has been compartmentalized by 
septum formation. Then, a long filament is produced, which consists of a chain of cells linked 
end to end [15]. Such a filament writhes as it elongates on the agar plate, and eventually 
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Figure 1: Snapshot of the filament configuration of bacterial cells on an agar plate at time 
t — 40 min after the first twofold part appeared. The scale bar indicates 20 /an. 


some segment starts folding and a twofold part of the filament is created. Figure [T] shows a 
typical configuration of a filament in which twofold parts have been created. Sooner or later, 
we will see the appearance of threefold parts, fourfold parts, and so forth, and the filament 
configuration of linked cells becomes complicated on the agar plate. 

Mendelson and coworkers have intensively studied the supercoiling processes performed 
by such cell filaments of B. subtilis They are interested in the situation 

wherein the filaments twist to make a double-stranded helix. The double-stranded structure 
itself twists while writhing, eventually comes in contact with itself, and forms a supercoil. 
By the repetition of such supercoiling processes, macroscopic structures of millimeter length 
are created, which are called bacterial macrofibers. The interesting motion of macrofibers 
was reported by Mendelson and coworkers. [23, 12T11231 [2b] , Kumada et al. also observed the 
growth of filaments of cells without separation for B. subtilis , in which the dependence of 
the morphology on C a was systematically studied [27]. Such growing filamentous cells were 
also observed for Escherichia ( E.) coli [28] . 

In this paper, we study the simplest situation wherein a filament of bacterial cells does not 
twist, and hence helical structures are not formulated. In our case, a single self-elongating 
filament repeatedly folds upon itself and shows a crossover from a one-dimensional structure 
to a two-dimensional structure. In Sect. [21 we explain the experimental setup and proce¬ 
dure for recording microscopic snapshots of cell filament configurations for about 6 h. We 
performed fractal analysis of filament configurations by the box-counting method. In Sect. 
[3l we report that the results imply fractal structures in the filament configurations, and the 
evaluated fractal dimension D shows a crossover from D = 1 to 2 as the growth process goes 
from the exponential phase to the stationary phase. Detailed study of the time evolution 
of the filament configuration of bacterial cells is discussed in Sects. [4] and [5j In Sect. [4J 
we propose a method of characterizing filament configurations on an agar plate by a set 
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of lengths of the simple part and the fc-fold parts with k = 2,3,4,, and the results of 
experimental measurements of these lengths are shown. To analyze the data given in Sect. 
SI we introduce systems of differential equations in Sect. 0and the nonlinear fitting of data 
to their solutions is performed. Section | 6 ] is devoted to concluding remarks. 

2 Experimental Procedures 

We observe a multiple-fission process from a single bacterial cell in the resting period under 
hard-agar and nutrient-rich conditions. The experimental setup and procedure are as follows. 

A solution containing 5 g of sodium chloride (NaCl), 5 g of dipotassium hydrogen- 
phosphate (K 2 HPO 4 ) and 10 g of Bacto-Peptone (Becton, Dickinson and Company, Franklin 
Lakes, NJ, USA) in 1 L of distilled water is prepared. The environmental parameter C n is 
given as the concentration of Bacto-Peptone; C n = 10 g-L _1 . Then, the solution is adjusted 
to pH 7.1 by adding 6 N hydrochloric acid (HC1). Moreover, the solution is mixed with 
10 g of Bacto-Agar (Becton, Dickinson and Company), which determines the softness of 
a semisolid agar plate. The environmental parameter C' a is given as the concentration of 
Bacto-Agar; C a = 10 g-L -1 . The environmental condition realized by these values of C a and 
C n gives a typical Eden-like pattern of B. subtilis colonies. The mixture is autoclaved at 121 
°C for 15 min, and 20 mL of the solution is poured into each sterilized plastic petri dish of 
88 mm inner diameter. The thickness of the semisolid agar plates is about 3.2 mm. After 
solidification at room temperature for 60 min, the semisolid agar plates are dried at 50 °C 
for 90 min. 

3 /iL of the bacterial suspension is inoculated at the center of each agar plate surface. 1 
fiL of the suspension includes about 10 2 cells in the spore state. The agar plates are left at 
room temperature for about 60 min to dry the bacterial suspension droplet. Thereafter, they 
are cultivated in a stage top incubator at 35 °C (INULG2-OTOR-CV, Tokai Hit, Shizuoka), 
which is attached to the stage of an optical microscope (1X71, Olympus, Tokyo). The spores 
at the inoculated center spot germinate about 2 h after the inoculation. They repeatedly 
undergo cell multiplications without splitting or movement, and form a long filament. The 
filament of cells linked in tandem grows two-dimensionally on the agar plate surface like 
a self-elongating string. This growth process of a filament of cells is observed through an 
optical microscope with a 20 x objective. 

A digital camera (DP71, Olympus, Tokyo) is connected to the optical microscope. Mi¬ 
croscopy snapshots are recorded every minute for about 6 h using bio-imaging analysis 
software (Lumina Vision, Mitani, Fukui and Tokyo) from the time when spore germination 
occurred to the time when the entanglement of a filament started two-dimensional colony 
expansion. 
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3 Crossover from One-Dimensional Structure to Two- 
Dimensional Structure 

In this section, the time after the inoculated spot germinated is denoted by T [min]. Let 
L[T) [pm] be the total length of a filament of bacterial cells observed at time T. If the 
specific growth rate is denoted by p, it will show an exponential elongation, 

L(T) = L(0)e /iT , T > 0. (3.1) 

First, we confirmed that the observed L(T ) data up to T = 200 min are well described by 
the single exponential function Eq. (13.Ill , and the specific growth rate was evaluated to be 

p = 3.58 x 10” 2 min” 1 . (3.2) 

This gives a cycle time (doubling time) of 

In 2 

t = -= 19.4 min. (3.3) 

p 

The snapshots of filament configurations at T = 125 and 245 min are shown on the 
left-hand side in Fig. [2] with a resolution of 2040 x 1536 pixels. They are analyzed by the 
box-counting method to evaluate the fractal dimensions of the filament configurations of 
bacterial cells on the agar plate. The procedure is as follows. At each time T, the snapshot 
picture is divided by squares (two-dimensional boxes) of linear size e, and then is counted 
the number Nt(e) of squares containing the pixels occupied by the filament of bacterial cell. 
If the configuration has a fractal structure, Nt(e) is scaled as 

N t (e) ~ e~ D(T) (3.4) 

with the fractal dimension D(T). We changed the value of e from 1 to 1536 pixels (from 
2.1 x ICE 1 to 3.2 x 10 2 pm in the real scale). As shown in the log-log plots given on the 
right-hand side in Fig. [21 the data for T = 125 min show a power law in the range of 3.5 

pm < e < 5.6 x 10 1 pm in the real scale, and those for T = 245 min show a power law in the 

range of 2.5 pm < £ < 2.1 x 10 1 pm in the real scale. The slopes —1.16 and —1.96 in these 
log-log plots in Fig. [2] provide the fractal dimensions _D(125) = 1.16 and .0(245) = 1.96, 
respectively. The result implies that the filament configurations have fractal structures. We 
evaluated the fractal dimensions of the filament at T = 45, 85,165, 205, 285, 325, and 365 
min and the results are plotted in Fig. [3] We found that the data can be fitted to the 
following sigmoid function: 

D(T) = ci tanh [a(T - T 0 )] + c 2 (3.5) 

with T 0 = 1.82x 10 2 min, a = 1.61 x 10” 2 min” 1 , ci = 5.1 x 10” 1 , and c 2 = 1.54. The 
time evolution of the multiple fission of bacterial cells from the exponential phase to the 
stationary phase is well described by the time dependence of the fractal dimension of the 
filament configuration on the agar plate. 
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Figure 2: (Color online) Upper left (lower left) picture: snapshot of the cell filament at time 
= 125 min (T = 245 min). The upper right (lower right) graph shows the log-log plots of 
r(e) versus e. We find linear regions in the log-log plots, and the fractal dimensions of the 
filament configurations are evaluated as .0(125) = 1.16 and .0(245) = 1.96, respectively. 
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Figure 3: Time dependence of fractal dimension D(T) of the filamentous cell configuration 
on the agar plate. The data are well described by the sigmoid function given by Eq. (13.51) . 

4 Folding Processes and Experimental Measurements 

In this section, we denote the observation time by t [min] instead of T [min], since we will 
set t — 0 at the time when the first folding of a filament occurs as explained below. 

Let L(t) [/mi] be the total length of a filament of bacterial cells at time t. Here, we write 
the specific growth rate as a, and then L(t) should obey the following differential equation: 

j t L(t) = aL(t). (4.1) 

On the semisolid agar plate, the linear growth is unstable, and a filament of cells starts 
folding after a short time. It repeats the folding processes and the configuration of the 
filament on the plate becomes complicated very rapidly. Over time, the multifold structure 
becomes dense and its region spreads over the plates. In this way, the configuration of a 
filament of cells shows a crossover from the one-dimensional structure to the two-dimensional 
structure, as shown in the previous section. 

In an interval of a filament, if bacterial cells are linked in tandem and form a curved 
line without folding, the interval is called a simple segment. At each time t, we consider 
the union of all simple segments in a filament and call it the simple part. We represent its 
length by l\(t) [/tm]. As illustrated by [1] in Fig. |4j we see the following process: 

Elementary Process [1] : elongation of a simple segment, 

and l\ (t) rapidly increases with time t. 

We set t — 0 at the time when the first folding of a filament is observed. In this 
experiment, t = T — 60 [min]. For t > 0, the following process occurs: 

Elementary Process [2] : creation of a twofold segment by folding of 

a simple segment, 
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Figure 4: Illustrations of elementary processes found in the self-elongating process of a 
filament of bacterial cells in the early stage. [1] Elongation of simple segment. [2] Creation 
of a twofold segment by folding a simple segment. Twice the length r 2 is added to the length 
l 2 . [3] Creation of a fourfold segment by elongation of a twofold segment. Four times the 
length r 4 is added to the length Z 4 . [4] Creation of a threefold segment by folding of a twofold 
segment on a simple segment. Three times the length r 3 is added to the length Z 3 . 

as illustrated by [2] in Fig. |4j The twofold part is defined by the union of all twofold segments 
in the filament, whose length is denoted by l 2 (t) [/iin]. By definition, Z 4 (t) = L(t),l 2 (t) = 0 
for t < 0, while li(t) — L(t) — Z 2 (Z), hif) > 0 for t > 0. Sooner or later, we will see threefold 
segments, fourfold segments, and so forth. We call the union of all A;-fold segments the 
fc-fold part, and write the length of the /e-fold part as 4(f) [/ini], k = 2,3,4,.... See the 
elementary processes [3] and [4] in Fig. U which create a fourfold segment and a threefold 
segment, respectively. We will attempt to characterize the time evolution of the filament 
configuration of cells using the set of lengths (h(t), Z 2 (t), h(t), ■ ■ ■ ) developing in time t. Since 
Xlfc>i 4(0 = T( 0 , the data are regarded as a time-dependent ‘partition’ of an exponentially 
growing length L(t). 

We found a critical time i* such that when 0 < t < t*, the whole filament of cells 
consists of only simple segments and twofold segments, while when t > i*, we observe the 
appearance of threefold and fourfold segments in a filament of cells and the configuration 
starts to become complicated. In this experiment, we observed 

f* = 45 min. (4.2) 

It corresponds to the time T* = t* + 60 = 105 min after the inoculated spot germinated. 
Note that it gives the time when the fractal dimensions D(T) of the filament starts to show 
its rapid increase in Fig. EJ 

Figure [5] shows a snapshot of a configuration of bacterial cells at t = 60 min. We obtained 
an enlarged photocopy of the snapshot picture and traced the filament segments by hand. 
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Figure 5: Snapshot of a configuration of the filament of bacterial cells at time t = 60 min. 
Each interval with a number k represents a fc-fold segment of the filament, where k = 2,3,4, 
and 5. The intervals without a number are the simple segments. The scale bar indicates 20 
pm. 

Each interval with a number k represents a fc-fold segment of the filament, where k = 2, 3,4, 
and 5. The intervals without a number are the simple segments. We used an opisometer, 
which is an instrument for measuring the length of arbitrary curved lines on a sheet. If the 
length of a segment indexed k is r*,, the length lkit) of the /e-fold part of the filament should 
be the sum of krk over all fc-fold segments. We measured (lk(t))k >l up to time t = 90 min. 
The results are listed in Table [U where the lengths of parts with k > 3 are summed and 
the values of h + (t) = Y2k>3^(t) are given. In the experiment, we fixed the held of vision of 
the optical microscope. Just before t — 70 min, a tip of the simple part ran out of the held 
of vision, and after t — 70 min, twofold segments may have been created out of our held of 
vision. Hence, the values of L and l\ at t — 70, 80, and 90 and those of I 2 at t = 80 and 90 
are not given in Table [H 

5 Analysis by Systems of Differential Equations 

5.1 Exponential growth of total length 

Equation (14.ip is solved as 

Lit) = Ae at . (5.1) 
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Table 1: Experimental data. 


t [min] 

L[/im] 

/i [/irn] 

h [n m] 

CO 

+ 

,3, 

0 

151 

151 

0 

0 

10 

228 

189 

39 

0 

20 

315 

247 

68 

0 

30 

485 

340 

145 

0 

40 

733 

442 

291 

0 

50 

1070 

620 

362 

90 

60 

1590 

810 

584 

197 

70 

— 

— 

944 

535 

80 

— 

— 

— 

1260 

90 

— 

— 

— 

2710 


Here, A = L(0) is the total length of the filament at t — 0 when the first folding occurs. 

As shown in Fig. [ 6 j the time dependence of L{t) is described by Eq. (15.ip very well for 
0 < t < 60 min, and we obtained the following values by the semilog fitting of the data: 

a = 3.95 x 1(T 2 min” 1 , A = 1.50 x 10 2 /jm. (5.2) 

The evaluation of a is consistent with the evaluation given by Eq. (13. 2 p for the specihc 
growth rate /i averaged over the longer time period 0 < T < 200 min. 

5.2 Systems of differential equations and their solutions 

In the time interval 0 < t < t*, only elementary processes [1] and [2] take place. We assume 
that in the elongation process of a cell filament with specific growth rate a, the ratio of the 
frequency of elementary process [2] to that of elementary process [1] is given by /3/(l — (3) 
with a constant 0 < /3 < 1. Then, the time evolution of the lengths li(t) and l 2 if) will be 
described by the following system of linear differential equations: 

^h(t) = “(1 ~P)h(t), 

^/ 2 (t) = al 2 (t) + a/3h(t), 0 < t < f*. (5.3) 

Note that the first term in the second equation in Eq. (15.31) describes the self-elongation 
process of the twofold part. Since we have set t = 0 at the time when the first folding occurs, 
Eq. (15.31) should be solved under the conditions 

/i (0) = L(0) = A, l 2 ( 0) = 0. (5.4) 

The solution is then given by 

h(t) = Ae a(1 -W, 

l 2 {t) = Ae a \ l-e” a/3t ), 0 <t<U. (5.5) 
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0 20 40 60 

Figure 6: Values of In L(t) plotted for t — 0,10,20,..., 60 min. The linear fitting of Eq. 
(15.11) in this semilog plot determines the values of a = 3.95 x 10” 2 min -1 and A = 1.50 x 10 2 
fi m. 

For t > t*, we take into account the following elementary processes in addition to pro¬ 
cesses [1] and [2]: 

Elementary Process [3] : creation of a fourfold segment by elongation of 

a twofold segment, 

Elementary Process [4] : creation of a threefold segment by folding of 

a twofold segment on a simple segment. 

As illustrated by [3] in Fig. SI a part of the elongating twofold segments becomes a fourfold 
segment. We assume that the ratio of the frequency of fourfold segment creation to that of 
simple elongation of the twofold part is given by 7/(1 — 7) with a constant 0 < 7 < 1. As 
illustrated by [4] in Fig. [4] the creation of threefold segments can occur only if a twofold 
segment touches a simple segment and folds on it and if they merge into a threefold segment. 
Thus, its frequency will be proportional to the product of li{t) and hit). We assume that 
this process reduces the total length of the simple part by <Mi(0^ (t) and that of the twofold 
part by S 2 h(t)l 2 (t) with transition rates per unit length Si > 0 and 82 > 0. Then, if we set 
I 3 + (t) = ^2 k>3 h(t) = h(t) + U(t) + ■ ■ ■ , we will obtain the following system of nonlinear 
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differential equations: 

= a;(l — /3)Zi(t) — <5iZi(t)Z 2 (t), 

= a(l - 7)4(0 + a/3h(t) - 8 2 h{t)l 2 (t), 

^4+(0 — a h+(t) + cr+(0 + (^i + <4)4 (t)l 2 (t), t > f*. (5-6) 

We assume that the parameters <4 and <4 are sufficiently small and solve the system of 
nonlinear differential equations given by Eq. (I5.6j) by perturbation. This assumption will be 
verified by the data fitting as explained in Sect. 15.31 

For k — 1,2, and 3+, we expand 4(f) as power series of <4 and 82 as 

4(0 = 


where l[ mi ' ,n2 \t) are time-dependent coefficients of the expansion. The TVth-order approxi¬ 
mate solution, Af = 0 , 1 , 2 ,..., is given by 

N 

C(‘) = E E s 7 'sri { r ,mi) (t)’ * = 1,2.3+. (5.8) 

72—0 7721 > 0,7722 > 0 , 

7721+7722=72 

In the following, we calculate the 0th- and first-order approximate solutions: 

it\t) = ~r\t\ 

l ( k\t) = l { k\t) + 8 l l^°\t) + 5 2 l^ 1 \t), k = 1,2,3+, t > f*. (5.9) 

The Oth-order approximate solution, (4° (f) : k = 1,2,3+}, solves the system of linear 
differential equations obtained from Eq. (15. 6 p by setting <4 — 82 — 0. In addition to the 
initial conditions corresponding to Eq. (I5.4[) . 

4 0 ) (0) = + 4 0 ) (0) = 0, (5.10) 

the definition of the critical time t* gives 

4">((.)=0. (5.11) 

Under these conditions, we have the following: 

4 0 ) (f) = Ae a(1_/,)t , 

4 0 ) (f) = -^_ e Q( 1 - 7 )t (l-e-“ (/ 3 " 7 ) 4, 

4° + } (f) = J^ eat {-/ 3 (e-“ 7 t -e-“ 7 '‘)+ 7 (e- Q/ 3 t -e- Q/3t *)}, t > 4 . (5.12) 


OO OO 


^2 s? i spit 1,ma) (t) 


7721 =0 7722=0 
OO 


7721 r7722 T( m l , m 2) 


E E ■W* 

72—0 7721 > 0,7722 > 0 , 

7721+7722=72 


(t), 


(5.7) 
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To express the first-order approximate solution, we introduce the multiple integrals 


J (1) (t;f*,a 1) = 

[ dse-° is 4 0) (s)4 0) (s), 

Jt , 


I {2 \t;U,a 2 ,a 1 ) = 

f dse _a2S / (1) (s; t*, a 3 ), 

Ju 


J (3) (f; £*, a 3 , a 2 , ai) = 

f dse _ “ 3S / (2) (s; t*, a 2 , ai), 

Ju 

(5.13) 


where a*, i = 1,2, and 3 are constants. By inserting lf ] (t) and l^\t) given in Eq. (15.12)) . 
they are calculated as 


i) 

J ( 2 ) (t;t*,a 2 ,ai) 


?(p—ai)t _ p(p—ai)t* 


3 (g-ai)t _ p{q—a-i)U 


f3A 2 

/? — 7 [ p — ai q — a i 

(3A 2 j'g(p-“l"“2) i _ g(p-“l-“2)(» g(9—ai-a2)t _ g(9—ai—Q2)*« 


/3 — 7 l (p - ai)(p - ai - a 2 ) 


(g - ai)(g - ai - a 2 ) 


1 


3 


3 (p-ai) 4 , 


•,(q-ai) 4 » 


/ ( 3 ) (t;t*,a 3 ,a 2 ,ai) = 


a 2 V p - ai 

®2 ^3 


q-a ± 

g(p a 2 ^3 )^* 


PA 2 

P~7 LI (p - <*i)(p- ai - a 2 )(p- ai - a 2 - a 3 ) 

g(<j—01—02—03)4 _ g(g—01—02—03)4* 

(<? - ai)(g - ai - a 2 )(g - a 3 - a 2 - a 3 ) 

g —(02+03)4 _ g —(02+03)4* /g(p—oi) 4 * g(g—ai) 4 * 


a 2 (a 2 + a 3 ) 


g 034 _ g 034 * 


3 —024* 


P - cq 

g(p—01—02)4* 


q ~ a\ 


g(cjf oi 02 ) 4 * 


a 2 


(p-ai)(p-ai-a 2 ) (<? - ai)(g - a 3 - a 2 ) 

g(p-oi) 4 * g(?-oi) 4 : 


P ~ «1 


q ~ ai 


(5.14) 


where p = a (2 — ft — 7 ), q = 2 a(l — /3). 

Then, the first-order approximate solution is given by 

4 1 } (f) = Zi 0 ) (t) — t*, a(l - ^)), 

4 1] (t) = 4°)(t) _5 ia/ 3e Q(1 - 7)i / (2) (t;t*,-a(7-^),a(l-/3)) 
-h 2 e“ ( 1 ~ 7 )t / ( 1 ) (t;t*,ai(l - 7 )), 

4?(f) = 4°>) 

+ 5i| - a 2 /hye at / (3) (£; £*, 07 , - 0(7 - /3), a(l - /?)) + e Qf / ( 1 ) (t; t*, a) j 
+ h 2 | - a 7 e at / ( 2 ) (t; t*, a, a(l - 7 )) + e" 4 / ( 1 ) (t; t*, a)|, t>t*. (5.15) 


Note that, owing to the initial condition 4^(0) = 4°4o) = A / ( /^(t) given in the first line of 
Eq. (15.15)) becomes independent of <f 2 . 
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Figure 7: (Color online) Data of h(t), l 2 (t), and Z 3+ (f) plotted by 0> A, and <), respectively. 
The thin curves represent the solution given by Eq. (15.5|) for 0 < t < i* with the parameters 
given by Eqs. (15.21) and (15.161) . The solution given by Eq. (15.151) for t > t* is represented by 
the thick curves, where the parameters are given by Eq. (15.171) . The fitting is excellent. 

5.3 Nonlinear fitting 

First, we used five pairs of data (h(t), hit)) of t — 0,10,..., 40 < t* = 45 in Table HI Note 
that the parameters a and A have already been determined using Eq. (5.2). By fitting to the 
solution given by Eq. (15. 5j) of the system given by Eq. (15.3j) of linear differential equations 
for (h(t), l 2 (t)), 0 < t < t*, we obtained the value of parameter 

13 = 0.317. (5.16) 

Next, we used the data 4(f), k = 1,2, and 3 + for t = 50, 60,..., 90 > t* = 45 in Table 
U] and performed their nonlinear fitting to the first-order approximate solution given by Eq. 
(I5.15|) of the system of nonlinear differential equations given by Eq. (I5.6|) . Here, A, a, and 
(3 are fixed to be the values given by Eqs. (15. 2 p and (15. 16j) . and 7, A, and 5 2 are chosen as 
fitting parameters. They are evaluated as 

7 = 0.313, 

<4 = 6.24 x 10“' min -1 • /inf 1 , S 2 = 3.75 x 10“ 8 min -1 • / um“ 1 . (5.17) 

The fitting is excellent as shown by Fig. [3 The evaluation given by Eq. (15 .1 7p is consistent 
with the assumption |<5j| <C 1, i = 1, 2, on which we solved the system of nonlinear differential 
equations given by Eq. (15.6j) by perturbation in Sect. 15.21 

The nonlinearity is very small but necessary in the fitting. To demonstrate it, we show 
the Oth-order approximate solution given by Eq. (15. 12j) as curves in Fig. |8j Here, we used 
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Figure 8: (Color online) Data of h(t), l 2 (t), and Z 3 + (t) plotted by 0> A, and <}, respectively. 
The curves show the Oth-order solution given by Eq. (15.121) . which ignores the nonlinear 
terms in Eq. (]5.6(). We failed to fit h+(t). 


the same values of A, a, and f3 as in Fig. [7J but we put Si — S 2 — 0. Figure [8] shows that if 
we ignore the nonlinear terms in Eq. (15. 6p . we fail to fit the data of Z 3 + (f),f > t*. 


6 Concluding Remarks 

In this paper, we have reported that the growth process of cells of B. subtilis under hard- 
agar and nutrient-rich conditions allows the realization of the dynamics of a self-elongating 
filament with sequential folding on a plane. Such a multiple-fission process without cell 
separation is commonly observed in the early stage of the growth process even if the agar 
concentration is changed, while the structure and motion of an entangled filament of cells 
in the later stage depend on environmental conditions [18] [21:, [3, IT] I2T] . Takeuchi et al. 
reported a study on filamentous cells of E.coli [28] ■ Environmental conditions to realize 
such self-elongation of cell filaments should be clarified by a systematic study of the early 
stage of bacterial growth processes. The classification of the morphology and dynamics of 
a long filament of cells depending on environmental conditions, time periods, spatial and 
geometrical restrictions, and so forth will be an interesting future problem. 

Here, we have focused on the simplest situation wherein a filament of cells simply repeats 
folding processes as it elongates and isotropically spreads over a two-dimensional plate. 
Note that Mendelson and coworkers have very extensively studied filamentous cell growth in 
the situation wherein supercoiling processes create helical macrofibers and their chiral self¬ 
propulsion motion is observed [IE] [19] [20] EU [22] [23] EU [25] [26]. Even in our simple case, it 
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seems to be highly nontrivial to provide a proper description of the filament configuration, 
which rapidly becomes complicated as it elongates with sequential folding embedded in a 
plane. In this work, we have proposed describing the global development by the time- 
dependent fractal dimension D(T) and the local folding processes by the time evolution of 
partitions (lk(t))k> 1 of the exponentially growing total length L(t) of the filament of cells, 
where k = 1 for the simple part and k > 2 for the fc-fold parts. 

The analysis discussed in Sect. 03 could be regarded as a mean-field-type approximation 
in the following sense. Let us consider a magnetic spin system on a lattice. In the mean- 
field theory, to describe a phase transition, we consider only the magnetization as an order 
parameter, which is obtained by averaging over spin configurations. The magnetization per 
spin m(T, H) at temperature T in an external magnetic field H is calculated by approxi¬ 
mating the correlated many-spin system by a single-spin system in a mean field generated 
by the surrounding spins, which is assumed to be proportional to m(T,H). The propor¬ 
tionality coefficient can be called the effective coupling constant J e g. In this way, we obtain 
the self-consistency equation for m(T , H) with the parameter J e g in addition to the external 
parameters T and H. If we want to compare the experimental data of a magnetization 
process of some material with the mean-field theory, the parameter J e g should be evaluated 
by some additional experimental observation. In this analysis of the filament configuration 
with folding, we summed lk(t) over k > 3 to define h + (t). By this reduction of variables 
from the series (lk(t))k> i to the triplet (h(t), hit), h + (t)), we obtained the finite systems of 
coupled differential equations given by Eq. (15. 3 ft for 0 < t < i* and by Eq. (15.6ft for t > t*. 
They involve the parameters a, (3, 7 , <5i, and 5 2 . We have evaluated these parameters as well 
as t* by experimental observations. 

As shown by Eqs. (5.16) and (5.17), the evaluated f3 and 7 have almost the same value. 
We can verify that the quantities given by Eqs. (5.12) and (5.14) have finite values in the 
limit 7 —> /3, and then our first-order approximate solution given by Eq. (5.15) is also valid 
in the case where [3 — 7 . We should note that if f} = 1/3 (7 = 1/3), the ratio of the 
frequency of twofold-segment (fourfold-segment) creation to that of simple elongation of the 
single (twofold) part is given by f3/{l — /3) — 1/2 ( 7/(1 — 7) = 1/2). We will continue our 
study to answer the question whether this result, (3 ~ 7 ~ 1/3, is universal. 

To improve the description, we need to carry out further studies on growing elastic 
filaments. Theoretical investigations can be found in the literature (29 ft 1301 l3Tft 132]. We hope 
that this experimental evaluations of parameters controlling the folding processes, which are 
described as l\ —> l 2 , l 2 —> U, and l\ +1 2 —3 > l 3 using our variables, will be useful for testing 
the validity of a possible theoretical consideration in the future. 
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